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Abstract 

Coulomb implosion mechanism of the negatively charged ion acceleration in laser plasmas is 
proposed. When a cluster target is irradiated by an intense laser pulse and the Coulomb explosion 
of positively charged ions occurs, the negative ions are accelerated inward. The maximum energy 
of negative ions is several times lower than that of positive ions. The theoretical description and 
Particle-in-Cell simulation of the Coulomb implosion mechanism and the evidence of the negative 
ion acceleration in the experiments on the high intensity laser pulse interaction with the cluster 
targets are presented. 



Recent developments in the ultra-intense laser pulse technology resulted in the intensity 
at the level above 10^^ W/cm^ [IJ. The ultra-intense laser pulse interaction with matter 
opens such new research fields as the fast ignition of inertial thermonuclear fusion [2], the 
charged particle beam acceleration for medical applications and the development of 
compact sources of high energy electrons, ions and photons [H 13 [6] (see also the review 
article |7j and literature quoted in). 

The negative ion generation has attracted a great deal of attention due to various ap- 
plications as the beam injectors for magnetically confined fusion devices [H], as the heavy 
ion fusion drivers [H], and for utilizing them in large scale ion accelerators (TU]. Negatively 
charged ion beams are typically generated via the charge exchange process during the pos- 
itive ion propagation through the alkali metal vapor. However, since the charge exchange 
cross section rapidly decreases as the incident energy exceeds keV level, it has been an im- 
portant issue to find an effective way to generate the high energy negative ions. Laser plasma 
provides a source of the negative ions with the energy in keV~MeV range pT| \n\ [TB] . In 
the experiment presented in Ref. [TTj, the negative ions with the energy about MeV have 
been observed when the water droplets are irradiated by the ultra-short and ultra-intense 
laser pulse (40 fs, / ~ 10^^ W/cm^). The negative ions were generated inside the laser 
plasma and are accelerated together with the positive ions. As far as it concerns the posi- 
tive ion acceleration in the laser plasmas, there were invoked several mechanisms including 
the sheath acceleration [T3j, the Coulomb explosion [T3], the shock acceleration [15], the 
radiation pressure dominant acceleration [17J, the "after-burner" [TS], etc. On the other 
hand, the acceleration mechanism of negative ions, according to our knowledge, has not yet 
been studied so far. 

In the present Letter we propose the Coulomb implosion mechanism of the negative 
ion acceleration. When a cluster target is irradiated by an intense laser pulse and the 
ponderomotive pressure of the laser light blows away the electrons, the repelling force of 
an uncompensated electric charge of positive ions causes the cluster Coulomb explosion. 
In the case of a multi-species cluster with a relatively small number of the negative ions, 
where possible mechanisms of negative ion creation are discussed in Ref. |19]. the electric 
field formed by the positively charged component accelerates the negative ions inward. The 
negative ions leave the target, passing through the center or bouncing its vicinity. Below we 
formulate the theoretical description and present the results of the Particle- in- Cell simulation 
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of the Coulomb explosion/implosion of a multi-species cluster. The evidence of the negative 
ion acceleration in the experiments on the high intensity laser pulse interaction with the 
cluster targets is presented as well. 

At first, we consider the Coulomb explosion dynamics of a multi-species cluster. The 
cluster comprises the positive and negative ion components. Within the framework of hy- 
drodynamics description the positive and negative ion motion is governed by the continuity 
equations: dtn± + V ■ {n±v±) = 0, the Euler equations: dt^r^ + (v-tV) v-t = ±€±^1/ M±, and 
by the equation for the electric field: V ■ E = 47r '^{^e±n±) . Here E is the electric field, the 
subscript (±) stands for positive and negative ions, n±, \±, ±e±, and M± are their density, 
velocity, electric charge and mass. These equations admit a self-similar solution with homo- 
geneous deformation, for which the densities are homogeneous, n± = n±{t), and the veloci- 
ties and electric field are linear functions of the coordinates: Vi^± = rhik^±{t) {rrikj^i)'^ 
Ei = eij{t)xj. Here summation over repeated indices is assumed. The matrix mjfc-|- is the 
deformation matrix, rhik^± is its time derivative and {mkj^±)~^ is its inverse matrix (see Ref. 
|2Uj). Now we assume a spherical, cylindrical or planar symmetry of the flow, for which 
the deformation matrix has a diagonal form with equal diagonal elements, k±{t). In this 
case, the velocity of Lagrange element is given by v± = k±x^ and of the Euler element by 
v± = k±x/k±, and the density is n± = no,±/kj_, with no,+ and no- being the initial values 
of densities of the positive and negative ions; d equals to the dimension such that d = 1, 
d = 2 and d = 3 correspond to planar, cylindrical, and spherical geometry, respectively. 
Substituting these functions to the above written equations, we obtain for k±{t) a system of 
two nonlinear ordinary differential equations: 
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where we introduced 00+ = A/47rno,+e+/M+, a;_ = ^y4'^^no-e'l/M_, u;± = 
Ay47rno,-e+e_/M+, and u;=p = A/47™o.+e+e_/M_. 

If the density of negative ions is small enough, no_ ^ e+no,+/e_, then we assume it 
vanishes, uq^ = 0, at the first step of approximation with respect to a small parameter 
C-'^'O -/e+no,+. In this case for d = 3, Eq. ([T]) describes the Coulomb explosion of a positively 
charged cluster [2T]. For initial conditions, fc+(0) = 1 and dk^/dt{0) = 0, its solution reads 
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a/A;+(A;+ — 1) + log ^a/AJ+ + \/kj^ — = -\/2c<j+t. In the case of cylindrical geometry, 
d = 2, which may correspond to the Coulomb explosion of the self-focusing channel |22] . 
Eq. (1) yields Erfi ^^/log kj^ = y^2/nLJ^t. Here Erfi (2;) = erf {iz)/i is the complex error 
function which is a real function of its argument. In a planar geometry, = 1, we have 
k+ = l + {uj+tf/2. 

Then, assuming the dependence of k+ on time to be given by the above written ex- 
pressions, we obtain from Eq. (2) the equation (Pk^/dt^ + Vt^_{t)k^ = uP_/k^^ with 
n^{t) = uj^/{k+{t)Y/^ and with initial conditions A;_(0) = 1 and dk_/dt{Q) = 0. If 
uj_/Q_ <^ 1, the negative ions during the first step move towards the center from where 
they bounce and move outwards. The ion velocity can be estimated to depend on the initial 
position, as l^;-! = [J^ VL^{t')dt'^ x". In the opposite limit, when uj-/VL- ^ 1, we seek a 
solution in the form of slowly varying and fast oscillating components, k^{t) = K^{t) + k^{t), 
with = 0, where < ... > denotes a time averaging. The component corre- 

sponds to the equilibrium solution of Eq. ^ for frozen dependence on time of i.e. 



K_{t) = (uj^/uJzfz)'^^'^ k^{t). For oscillating component, k_, we find the equation of 
oscillator with a time dependent frequency: d^k^/dt^ + dVt^_{t)k_ = 0. Due to a small- 
ness of the ratio uj-/uj^ ^ 1 a dependence on time of VL^{t) is slow and we can use the 
WKB approximation. This yields the expression for the oscillating part of k^{t). It reads 
k^(t) = [A;+(t)]'^''^ cos Jp r2_(t')(it'j , which describes the particle oscillations with grow- 
ing amplitude and decreasing frequency. 

The ion energy spectrum can be found by considering a dependence of the velocity on 
the ion initial position. We see that the ion kinetic energy, S = M_f^/2, is proportional 
to the square of the ion initial coordinate: S oc (a;°)^. Using a condition of the continuity 
in the phase space, dN^ = A7mo^^{x^Y^^dS /\dS/dx^\ we find that the ion energy spectrum 
is inversely proportional to the square root of energy, dN/d£ oc in the case of planar 

geometry; it is flat, dN/dS = constant in cylindrical geometry, and is proportional to the 
square root of energy oc ^/S in the case of spherical geometry. The negative ion energy is in 
a factor k = (lu_/lu^)^^'^ = (e_no,_/e+r;,o,+)^''°' smaller than the energy of the positive ions. 

Now we demonstrate the implosion dynamics of negative ions by using the two- 
dimensional Particle-in-Cell (PIC) simulations. The simulation conditions are as follows. 
The target is a cluster with the diameter of 1 fim composed of electrons, protons and neg- 
ative hydrogen {H~) ions whose densities are LOuc, Line, and OAuc, respectively. Here 
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FIG. 1: Radial density profiles of (a) protons and (b) negative hydrogen, H , ions during the 
Coulomb explosion and implosion. The densities are normalized on the critical density. 



FIG. 2: (a) Distribution of the radial component of electric field Er at time t=130 fs. The electric 
field is normalized on the incident laser electric field, which is 27TV/m in our simulations. The 
profile of along the horizontal line passing y=20 is shown in the lower part, (b) Trajectories of 
H ions. 

ric = rUeUj'^ /Ane"^ is the critical density for the laser light with the frequency u. In this 
case the parameter k equals 0.1. The electron initial temperature is set to be equal to 500 
eV. The ions are assumed to be initially cold. The target is located at the center of the 
simulation box which has a size of 40/im in both the x— and directions. The laser pulse 
irradiates the target from the left hand side boundary. It propagates in the x— direction and 
is polarized along the y— direction. The pulse has a fiat top form rising up in a laser period 
and keeping its peak intensity of 1.0 x 10^*^ W/cm^ which is constant for the duration of 15 
fs, where laser wavelength is 1 fim. The laser pulse is focused to the spot with the diameter 
of 3yum (FWHM), where the focal point is a; = 18/im. The simulation box has 4000x4000 
meshes, with ~ 1 x 10^ quasi-particles used. 

By the irradiation of the ultra-intense laser pulse, the electrons are heated and swept 
away from the target leaving the target positively charged. This results in the Coulomb 
explosion of the positive ions and in the Coulomb implosion for the negative ions. The time 
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history of radial density profile of protons and H~ ions along x-axis (x<20) is plotted in Fig. 
[T] (a) and (b), respectively (r = — x+20). The protons expand radially due to the Coulomb 
explosion, having a rather flat distribution as in the self-similar solution of two-dimensional 
case. H~ ions start to move inward from the outer region, and they are most compressed at 
t = 130 fs whose center is located at (x,y) = (19.95,20.0), which is 0.05 /im left to the target 
center. This off-center of the implosion comes from the asymmetry of laser irradiation. As 
the H~ ions are accelerated towards the center, the electric field, which sign is opposite to 
that of the outer region, is induced at the center due to accumulation of the negative charges, 
as it is seen in Figj2] (a). The electric field decelerates the imploding H~ ions and attract 
protons towards the center which is seen in Figjl|^a). The trajectories of H~ ions are shown 
in Figj2]^b). Particles with numbers from 1 to 3 are initially located at the periphery of the 
target. They are accelerated towards the center and slightly deflected by the electric field 
and leave the target. They have maximum energy of 1 MeV when they reach the cluster 
center, and then they are gradually decelerated while leaving the target. Particle number 4 
moves towards the center and then is reflected back by the induced strong electric field. Its 
energy at the time when it has left the target is 0.1 MeV. Particle number 5 being located 
near the center bounces between the outer and inner regions and then leaves the target. It 
achieves the energy being only ~30 eV. As a result, energetic H~ ions are generated from 
the periphery of the target. 

The Coulomb explosion/implosion of positive/negative ions is illustrated in Fig. [3| where 
the phase planes {y,Py) are plotted. In the row (a,b) we show the ion phase planes at time 
t=80 fs. We see the linear (self-similar) dependence of the ion momentum on the coordinate 
corresponding to the positive ion expansion and to the negative ion collapse. In the row 
(c,d) we plot the phase planes at the time t=180 fs, after the negative ions have bounced 
from the central region. We see that both the positive and negative ions expand radially. 
A local compression of the positive ions induced by the electric field formed at the center is 
also seen in Fig. |3|:;. 

The energy spectra of protons and H~ ions are shown in FigsQa) and (b), respectively. 
The spectra of protons are flat, which is in agreement with the theoretical prediction in the 
case of two-dimensional Coulomb explosion. The maximum energy of negative hydrogen 
H~ is equal to 1.0 MeV. We see that the negative ion energy is approximately ten times less 
than the energy of positive ions. This is in agreement with the above formulated theoretical 
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FIG. 3: Phase planes {y,Py) of the positive (in the frames (a) and (c)) and negative (in the frames 
(b) and (d)) ions. The rows (a,b) and (c,d) correspond to the time t=80fs and t=180fs, respectively. 
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FIG. 4: Time history of the energy spectra of (a) protons, and (b) H ions. 

model with = kS^, for the parameter k equal to 0.1 in our simulations. 

In the case of lower values of the laser intensity or/and higher values of total initial 
number of the ions, the laser light does not blow away all the electrons. For example, when 
a laser intensity is / = 1.0 x lO^^W/cm^ and a target has the same parameters as above 
calculations, we do not have a regime of pure Coulomb explosion. Nevertheless, the negative 
ions are accelerated towards the target center. Since the sheath field inside the expansion 
front is weakened by the bulk electrons, the energy achieved by negative ions does not 
oscillate and they sustain the energy after leaving the target. 

Negative ions are observed in experiments using water droplet pT] . or thin foil target [T2]. 
In our experiments negative ions are observed for the first time using mixtures of CO2 
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FIG. 5: A pattern registered on CR-39 with the Thomson parabola. 

clusters and He gas target. The laser has energy of 130 mJ with pulse duration of 35 fs. The 
pulse is focused into the spot with 30 /im in diameter, resulting in the intensity of 7 x 10^^ 
W/cm^ in vacuum. The CO2 clusters with the diameters equal to 0.35 /im, which contain 
5 X 10® molecules each, are generated using specially designed supersonic gas jet nozzle |23]. 
Average distance between the clusters is 5 /im. The laser pulse contrast of 10"^ at the ns 
time scale is expected to lead to the cluster heating and evaporation by the pre-pulse. This 
results in the interaction of the main, fs, laser pulse with the plasma clouds of larger volume 
and lower density than in the initial clusters. The registered parabolic line of the positive 
and negative ions accelerated in the Coulomb explosion/implosion are shown in Figjs} The 
Thomson parabola is positioned in the direction of 135° from the laser axis, i.e. we observe 
the ions moving in the backward-aside direction with respect to the laser beam propagation. 
In Figg we see the parabolic lines produced by negative C ions as well as the carbon ions 
up to C^"^ ion and by the oxygen ions up to O^"*" ions. The maximum energy of positive C^"*" 
ions is 4.8 MeV, and that of the negative C~ ions is 0.6 MeV. The ratio of maximum energy 
of positive and negative ions is 1/8, which corresponds to k = 1/8 with no,-/'n.o,+ — 0.02. 
The detailed discussion of experimental results will be presented in the forthcoming papers 

In conclusion, the proposed Coulomb implosion model explains the acceleration of neg- 
atively charged ions in the laser-plasma. The negative ions initially located at the target 
periphery are accelerated more efficiently. A final energy of the negative ions is several times 
less than the positive ion energy. The Coulomb implosion mechanism is clearly demon- 
strated by the PIC simulations. Acceleration of positive ions via the Coulomb explosion is 
a well-known mechanism. Negative ions can be accelerated in the same field, however, their 
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acceleration occurs in the opposite direction with the bouncing-back in the vicinity of the 
center of symmetry, which results in the generation of high energy negative ions. We note 
here that the mechanism is applicable in the case of thin foil, self-focusing channel, and 
cluster according to in ID, 2D, and 3D configuration. 
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and Technology commissioned by MEXT of Japan. The simulations are performed by using 
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